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1. Introduction 

As a result of the success of current Imaging Atmospheric Cherenkov Telescopes (lACTs) 
and the improvement of the different technologies involved, the next generation of ground-based 
Very High Energy (VHE) y-ray detectors is under development. The Cherenkov Telescope Array 
(CTA)^ [1] will scrutinize the y-ray sky from some tens of GeV up to more than a hundred TeV 
with unprecedented capabilities, improving the sensitivity of current instruments by up to an order 
of magnitude. 

In order to achieve these goals, the CTA Observatory will consist of two different sites, one in 
each hemisphere, and telescopes of three different size classes: Large Size Telescopes (LSTs) sensi¬ 
tive to the faint-low energy showers (below 200 GeV), Medium Size Telescopes (MSTs) increasing 
the effective area and the number of telescopes simultaneously observing each event {multiplicity) 
within the CTA core energy range (between 100 GeV and 10 TeV) and Small Size Telescopes 
(SSTs, for CTA South only) spread out over several km^ to catch the rare events at the highest en¬ 
ergies of the electromagnetic spectrum (up to ^ 300 TeV). The current design foresees 4 LSTs, 25 
Davies Cotton (DC) MSTs, and 70 SSTs for CTA South as well as 4 LSTs and 15 MSTs for CTA 
North. The southern array may be augmented with a proposed extension of up to 25 Schwarzschild- 
Couder (SC) MSTs. 

The proposed designs for the northern and southern observatories will provide full sky cov¬ 
erage with an improved sensitivity alongside better angular and energy reconstruction. The CTA 
southern site is expected to be larger, composed of ^ 100 telescopes, to take advantage of its priv¬ 
ileged location to observe the Galactic Plane and Galactic Center over the full VHE range. The 
northern site is expected to complement it with a reduced number of telescopes of the order of 
- 20 (LSTs and MSTs). 

In this work we present results from the second large-scale MC production (Prod-2), and 
provide a realistic estimation of the future observatory capabilities. Prod-2 results concerning the 
site selection (such as the performance dependencies on altitude or the local Geo-magnetic field) 
and the performance of mixed MST type layouts are presented as separate contributions [2, 3]. 

2. Monte Carlo simulations and analysis tools 

CTA Monte Carlo simulations, as detailed in [4], are performed by defining a large telescope 
layout, comprising few hundreds of telescopes of different types arrayed over an area of around 
6 km^. Simulations of extensive air showers (LAS) initiated by y-rays, cosmic-ray nuclei and 
electrons are generated using CORSIKA [5] together with the simulation of the telescope response 
using sim_telarray [6]. By selecting sub-samples of telescopes of approximately equal costs 
from the macro-layout and evaluating their performance, conclusions on the most cost-effective 
CTA array are inferred. 

The first large-scale MC production (Prod-1) already simulated a wide variety of conditions 
(using a constant geomagnetic field strength): two different altitudes (2000m and 3700m), obser¬ 
vations at zero, medium and large zenith angles (0°, 20° and 50°) and partial moon light. The main 
conclusions (as detailed in [4]) were the following: Several of the proposed layouts could satisfy 

^http://www.eta-observatory.org/ 
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CTA scientific expectations for most physics cases. In addition, the improvement at the lowest 
energies (E < 100 GeV) of high-altitude sites does not compensate the effective area loss above 
those energies, favouring sites at moderate altitudes. 

The second large-scale MC production (Prod-2) intended to go one step further, with clear 
objectives: to assess the impact of the construction site on the CTA performance and to search for 
a feasible layout candidate, taking into account updated cost estimates and new proposed telescope 
types. With this in mind, telescope responses, site properties and the layout composition were 
configured with more realistic parameters: 

Telescope response : Using all s i m_t e 1 a r r a y improvements and the updates from all the tele¬ 
scope designs (such as Photo-Multiplier Tubes efficiencies, mirror reflectivity or readout systems), 
full ray tracing was performed, storing long readout windows per pixel to allow testing of cus¬ 
tom trace integration schemes in the analysis and including 3 different trigger schemes: majority, 
analog-sum and digital-sum. 

Layout composition: New large macro-layouts of telescopes were defined (one for each simu¬ 
lated hemisphere) with new telescope positions and types: 229 telescope positions in the southern 
site with up to seven telescope types were simulated (LSTs, DC-MSTs, SC-MSTs, and up to four 
variants of SSTs). For northern sites three telescope types (LSTs, DC-MSTs and SSTs) were sim¬ 
ulated at a total of 61 positions. 

Site properties: For each candidate site, specific atmospheric density profiles, altitudes and geo¬ 
magnetic fields (direction and intensity) are used. A total of 3 Northern and 6 Southern Hemisphere 
sites were simulated. 

Simulated showers include gammas (both point source and diffuse) and background (protons, 
nuclei, and electrons), with protons (^ 100 billion events per site) being the particle type consuming 
most of the CPU time and disk space resources, even though few of them trigger and pass the 
gamma selection cuts. Most of the simulations are produced for zenith angles of 20°, with showers 
coming from North and from South to include the impact of the geomagnetic field. In addition, 
dedicated productions were devoted to increased levels of Night Sky Background or to different 
telescope separations. The Prod-2 simulations were carried out on the CTA computing Grid as 
well as on the MPIK and DESY computer clusters. 

3. Analysis 

To evaluate the performance of a telescope layout, analysis tools similar to the ones used in 
the current lACT experiments are employed. Results presented here use methods based on the 
MAGIC Analysis and Reconstruction Software (MARS) [7], the software package used by the 
MAGIC collaboration. The analysis performs the following steps: 

• Trace integration: Pixel charge and the time of arrival are extracted from sim_t el array 
output for each triggered event. This is performed in a 2-step process, in which a fit to the 
photon arrival times for the most significant pixels is subsequently used to search for dimmer 
signals in a reduced time window. 
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• Image cleaning and parameterization: A traditional Hillas-like parameterization of the im¬ 
ages is then performed. 

• Stereo reconstruction: Using the information gathered by all triggered telescopes observing 
an event, a stereoscopic reconstruction is performed that measures the energy and direction of 
the primary gamma ray. Multivariate event classification algorithms are used for background 
suppression, discerning between cosmic and y-ray initiated showers. 

• Performance estimation: Optimal cuts are determined in order to calculate the resulting per¬ 
formance, expressed by the Instrument Response Functions (IRFs). Using the shower’s ar¬ 
rival direction and telescope multiplicity, differential sensitivity is maximized for each en¬ 
ergy bin. Sensitivity is computed requiring five standard deviations (5 < 7 ) at each energy bin 
(equation 17 from [8] considering a ratio of background to signal exposure of 5). In addition, 
the excess needs to be larger than 10 y-rays, and must correspond to a rate higher than the 
5% of the residual background rate. 


It is worth to note that alternative analysis 
tools have been developed for the analysis 
of CTA MC data, based on different soft¬ 
ware packages belonging to different lACT 
experiments. Fig. 1 shows the differential 
sensitivity of the “2Q” layout in 50 hours 
calculated with 3 different analysis chains 
(all introduced in [4]): The Baseline analysis 
at the MPIK, the evndisplay based analysis 
at DESY and the MARS based one detailed 
here. Although some differences appear, re¬ 
sults are consistent taking into account the 
differences between analyses (image clean¬ 
ing, shower reconstruction, quality cuts and 
background rejection power). 



Figure 1: Differential sensitivity for the CTA-S 
“2Q” candidate array (50 hours of observation, 
N/S pointing average) calculated with 3 alterna¬ 
tive analysis chains: Black: MARS analysis at 
IFAE. Red: Baseline analysis at MPIK. Green: 
Evndisplay analysis at DESY. 


4. CTA performance 

This work focuses on describing the capabilities of two feasible CTA-N and CTA-S layouts. 
Here, the layouts considered will be the “2N” candidate array (4 LSTs and 15 MSTs) for the CTA- 
N layout simulated at Tenerife and the “2Q" (4 LSTs, 24 MSTs and 72 double mirror SC-SSTs) 
for the CTA-S simulated at Namibia (Aar), both shown in Fig. 2. The choice of sites is not related 
in any sense with the final construction location, and is based on availability of the simulated data 
files and the number of cross-checks with alternative analysis chains. 
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Figure 2: Proposed layouts for the Northern and Southern Hemisphere sites. Left: CTA-N candi¬ 
date layout “2N", made up of 4 LSTs and 15 MSTs. Right: CTA-S candidate layout “2Q” composed 
of 4 LSTs, 24 MSTs and 72 4m Schwarzschild Couder (SC) SSTs. 


Each CTA layout performance is calculated from the average of two simulated observations: 
north and south pointing directions at 20° of zenith angle. This is done to account for the per¬ 
formance differences coming from the geomagnetic field (GF) effect, detailed in [2]. It should be 
noted that the simulated observation directions for the chosen Southern Hemisphere site are partic¬ 
ularly optimistic, as the perpendicular component of the geomagnetic field for the north pointing 
direction is close to zero. Differences in the differential sensitivity caused by the GF between 
CTA-N and CTA-S decrease by ^ 35% when considering a realistic distribution of observation di¬ 
rections in zenith and azimuth, keeping in mind that Northern Hemisphere sites have, on average, 
higher perpendicular GF components. 

The minimum detectable flux in 50 hours from a steady y-ray point-like source as a function of 
the energy (differential sensitivity) for both CTA site candidates is shown in Fig. 3, together with 
the current MAGIC [9] and VERITAS^ (during 50 hours) and future HAWC [10] and Fermi-LAI 
[11] sensitivities. Both sites attain good performance over four decades of energy, and would be 
able to measure a spectrum 100 times fainter than the Crab Nebula (10 mCrab) across an order of 
magnitude in energy in 10 hours of observation for the case of CTA-S, increasing to 25 hours for 
CTA-N. 

As shown in Fig. 4, the CTA Observatory will not only provide an improvement in sensitivity 
with respect to the current instruments; it will also improve the angular and energy resolution, 
shrinking by a factor ^ 3 the radius of the Point Spread Function (PSF) with respect to current 
instruments, resolving structures on the order of 2 arc minutes at the highest energies observed 
by CTA, becoming the y-ray detector with the best angular resolution ever constructed. Energy 
resolution will also improve significantly, reaching values close to the 5% in fractional uncertainty. 

To understand the overall performance of such a system, it is required to recognize the con¬ 
tribution of each telescope type at different energy ranges. Performing independent analyses of 
sub-layouts of each telescope type within the CTA-S candidate array, it is possible to see the im¬ 
pact of each telescope on the differential sensitivity. Fig. 5 shows the detached performance of 

^http://veritas.sao.arizona.edu/ 
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Figure 3: Differential sensitivity for a 
point-like y-ray source of the CTA-N 
and CTA-S candidate arrays (50 hours 
of observation, N/S pointing average) 
together with the current MAGIC [9] 
and VERITAS (50 hours) and and fu¬ 
ture Fermi-hAI [12] (over 10 years of 
operation) and HAWC [10] (1 and 5 
years) attained sensitivities. 

10 '^ 10 '^ 1 10 10 ^ 

Energy (TeV) 





Figure 4: Angular and energy resolution for both CTA-N and CTA-S candidates with respect to the 
reconstructed energy: {Left) The angular resolution, defined as the angle containing the 68% of the 
reconstructed gamma-rays, relative to the true direction. (Right) The energy resolution is defined 
such that 68% of gamma rays will have true energy within AE of their reconstructed energy. 


each layout of individual telescope types contained in the CTA-S candidate: 

• Below 100 GeV: The subset of 4 LSTs dominates the sensitivity below 100 GeV. Although 
this range may seem small, a large number of physics cases will benefit from the CTA per¬ 
formance in this region, and cross-calibration with other y-ray detectors such as Fermi-LAI 
is desired. 

• From 100 Gev to 5 TeV: The layout of 24 MSTs dominates the instrument performance in 
the core energy range, above 100 GeV up to 5 TeV. Their contribution to the lower energies 
is also significant, as they allow to reconstruct events of ^ 60 GeV (boosting sensitivity by 
^ 20%). They also improve gamma-hadron separation of events observed by the LSTs, even 
at the lowest energies, as they detect possible ttq initiated electromagnetic sub-showers or 
muons originated in hadronic EASs. They also contribute to higher energies up to ^ 10 TeV. 

• Above 5 TeV: The array of 72 SC-SSTs, covering an area of ^ 4 km^, dominates the sensi- 
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Figure 5: Different telescope type con¬ 
tributions to the over-all differential 
sensitivity of CTA-S candidate layout 
“2Q", simulated at the Namibian site 
(50 hours, north/south pointing aver¬ 
age). LSTs govern the low energy 
range below 100 GeV. MSTs dominate 
in the CTA core energies, up to 5 TeV, 
where SSTs start to monopolize perfor¬ 
mance. 

- 1.5 -1 - 0.5 0 0.5 1 1.5 2 

log^^(Energy/TeV) 

tivity above 5 TeV (CTA-S only). Comparison of alternative layouts included in the Prod-2 
simulations shows that the performance of about 70 SC-SSTs (of 4 m diameter) is compara¬ 
ble to that of about 35 DC-SSTs of 7 m diameter when both arrays cover a similar area. 

Off-axis capabilities are also crucial for a significant number of the CTA key science projects. 
Sky surveys and studies of diffuse and extended sources of gamma-ray emission will benefit from 
a wide field of view (FoV) and accurate event reconstruction away from the camera center. To 
characterize the observatory performance for different off-axis angles, the differential sensitivity of 
a point-like source located at different distances from the center of the camera is estimated. 

Below 100 GeV, performance is dominated by the LSTs, which are telescopes with a smaller 
FoV. Their sensitivity does not decrease significantly at distances from the camera center smaller 
than 2°, but declines sharply at larger distances (by a factor 4 at a 2.5°). For the MSTs, their 
modified Davies-Cotton optics improve off-axis performance with respect to the parabolic one. 
Between 100 GeV to 5 TeV, the sensitivity drops by less than a factor 2 on sources 3° away of the 
center of the camera. Above 5 TeV the off-axis performance improves significantly. This effect 
is related to the good off-axis performance of the SC-SSTs dominating at these energies. Their 
performance is expected to be approximately flat up to 5° away of the camera. 

5. Conclusions 

The analysis presented in this work of both feasible CTA-N and CTA-S candidates represent a 
detailed estimation of the future CTA capabilities, capable of fulfilling all scientific goals defined 
by the Consortium. With the results from both CTA large-scale MC productions, layouts have 
been optimized to find the best trade-off between event quantity and reconstruction quality. The 
search for the most efficient layout for each site is still on-going, and will likely be completed in 
the third large-scale MC production (Prod-3), which is currently in preparation. A denser layout 
of 4m-SSTs is capable of attaining similar performance as the one initially proposed with 7m- 
SSTs. In addition, the newly designed Schwarzschild-Couder SSTs show great off-axis capabilities 
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providing, together with the MSTs, a FoV with a radius larger than 3° for the whole energy range 
above 100 GeV, increasing up to ^ 5° at the highest energies. 

MC simulations provide a crucial test-bench for the different designs within the CTA project, 
and Prod-2 results demonstrate that their correct implementation would attain the desired perfor¬ 
mance and potential scientific output. 

Acknowledgments 

This work is partially funded by the ERDF under the Spanish MINECO grant FPA2012-39502 and 
FPA2010-22056-C06-06. We gratefully acknowledge support from the agencies and organizations listed 
under Funding Agencies at this website: http://www.cta-observatory.org/. 

References 

[1] M. Actis and et al.. Design concepts for the Cherenkov Telescope Array CTA, Experimental 
Astronomy 32 (Dec., 2011) 193-316, [arXiv: 10 08.3703]. 

[2] G. Maier et al., Monte Carlo Performance Studies of Candidate Sites for the Cherenkov 
Telescope Array, . 

[3] T. Hassan et al.. Layout design studies for medium-sized telescopes within the Cherenkov 
Telescope Array, . 

[4] K. Bernlohr et al., Monte Carlo design studies for the Cherenkov Telescope Array, 
Astroparticle Physics 43 (Mar., 2013) 171-188, [arXiv: 1210.3503]. 

[5] D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz, and T. Thouw, CORSIKA: a Monte Carlo 
code to simulate extensive air showers. Forschungszentrum Karlsruhe GmbH, Feb., 1998. 

[6] K. Bernlohr, Simulation of imaging atmospheric Cherenkov telescopes with CORSIKA and 
sim_ telarray, Astroparticle Physics 30 (Oct., 2008) 149-158, [arXiv :0808.225 3]. 

[7] A. Moralejo et al., MARS, the MAGIC Analysis and Reconstruction Software, ArXiv e-prints 
(July, 2009) [arXiv:0907.0943]. 

[8] T.-P. Li and Y.-Q. Ma, Analysis methods for results in gamma-ray astronomy, Astrophysical 
Journal 111 (Sept., 1983) 317-324. 

[9] MAGIC Collaboration, The major upgrade of the MAGIC telescopes, Part II, ArXiv e-prints 
(Sept., 2014) [arXiv:1409.5594]. 

[10] A. U. Abeysekara et al.. Sensitivity of the high altitude water Cherenkov detector to sources 
of multi-TeVgamma rays,AP 50 (Dec., 2013) 26-32, [arXiv: 130 6.580 0]. 

[11] M. Ackermann et al.. The Fermi Large Area Telescope on Orbit, Astrophysical Journal 
Supplement 203 (Nov., 2012) 4, [arXiv: 120 6.18 96]. 

[12] S. Funk, J. A. Hinton, and CTA Consortium, Comparison of Fermi-IAT and CTA in the 
region between 10-100 GeV, AP 43 (Mar., 2013) 348-355, [arXiv: 1205.0832]. 


8 



